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Abstract 

We show how the generation of a lepton number in the Early Uni- 
verse induced by Va ^ i^sj in presence of smaU baryon number in- 
homogeneities, gives rise to the formation of lepton domains. The 
structure of these domains reflects the spectral features of the baryon 
number inhomogeneities that generated it and interestingly the size 
of the lepton domains can be super-horizon. 
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1 Introduction 

If one assumes the existence of a mixing between an active neutrino flavor a 
and a sterile neutrino flavor, then a lepton number can be generated in the 
early universe [0 . To this purpose a small vacuum mixing angle (sin 26q <C 
1) is necessary (the mass eigenstates almost coincide with the interaction 
eigenstates). Moreover the sterile neutrino must be lighter than the active 
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neutrino and the absolute value of (5m^, the difference of squared masses, 
must be in the range (10~^ -7- 10^) eV^. In fact the critical temperature at 
which the generation occurs is given, if the sterile neutrino production is until 
that time negligible, by the expression: 



It follows that if \Sm'^\ ~ 10 ^eV^, then Tc ~ 3 MeV and the lepton number 
variation would be dominated by the MSW effect. In this lepton 
number larger than 10"'' cannot be generated 0, |^. 

On the other hand if |5m^| ~ lO^eV^ then Tc ~ ISOMeV and one should 
be able to describe the neutrino oscillations in a quark-gluon plasma, some- 
thing beyond the present level of matter effects account (moreover it is not 
interesting in this context, because it would imply active neutrino masses 
much higher than 100 eV, cosmologically excluded if neutrinos are stable). 

A rigorous description of neutrino oscillations in the early universe re- 
quires a quantum kinetic approach able to describe the evolution of the sta- 
tistical density matrix for the two mixed neutrino flavors . Such a descrip- 
tion must include the different effects of matter on the neutrino mixing: a 
coherent effect due to forward scattering a loss of coherence due to the 
collisions that change the neutrino momentum [^, |^ and the repopulation of 
active neutrino states (depleted by the oscillations into the sterile neutrino 
states) through the collisions ||, |^ . 

For a better understanding of the generation of a lepton number, it has 
been shown to be more convenient to turn to a simpler description that ne- 
glects the possibility of a MSW effect at the resonance. This approximation 
proves to be valid at temperatures T ~ 5 MeV because the collisions de- 
stroy the coherence at the resonance. In this way the physics underlying the 
lepton number generation is isolated and can be clearly understood PH] ]. 
The role of collisions, from this point of view, is crucial: they can be use- 
fully considered, in a heuristic sense, as a measurement process able to make 
the two-quantum state collapsing in one of the two interaction eigenstates. 
Through this effect sterile neutrinos are produced and the presence of a small 
baryon number induces a tiny asymmetry between the production rates of 
neutrinos and antineutrinos. If the sterile is lighter than the active neutrino 
and if the values of mixing parameters are in the intervals previously indi- 
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cated, this initial asymmetry is amplified through the generation of a lepton 
number that starts at the critical temperature Tc. 

In this mechanism the momentum dependence plays an important role 



11] and a monochromatic approximation provides only a rough descrip- 



tion. The interesting feature is that, within this simplified physical picture, 
the quantum kinetic description through the density matrix collapses into 
an effective (classic) kinetic description (Pauli-Boltzmann approach). This 
is possible because in the regime where collisions are rapid enough, the evo- 
lution of the off-diagonal terms can be disentangled by the diagonal terms 
evolution and explicit expressions can be derived for them. In this way, in the 
description of the evolution of the statistical properties of mixed neutrinos, 
only the diagonal terms are left, the usual statistical distributions: a simpler 
physical picture yields a simpler mathematical description. The derivation 
of the equations for the distributions can be done either directly from the 
simplified physical model |T0[ or also, with a more formal procedure, from the 
quantum kinetic equations themselves, via some approximations, indicated 
by the authors globally as static approximation, valid under appropriate con- 



ditions on the mixing parameters [12|. The differences that one expects from 



a full quantum kinetic description are at low temperatures (T ~ 5MeV), 
when the MSW effect becomes important and must be taken into account, 
and at the critical temperature for large enough vacuum mixing angles, when 
the growth of lepton number is so rapid that an adiabatic condition to de- 
scribe the process does not hold anymore. This adiabatic condition is the 
possibility to neglect any change in the effective potentials, and thus in the 
mixing angle, between the collisions on average. The collisions, if vacuum 
mixing angles is not too large, are able to average out the coherent effects of 
the oscillations in the macroscopic quantity (like the neutrino asymmetry). 

This analysis on the validity of the approach is confirmed by the numer- 
ical calculations performed using the QKE. In [|1^, |14| it has been shown 
how the MSW effect at low temperature is able to amplify the growth of 
lepton number up to values slightly lower than the maximum absolute value 
obtainable of 3/8, corresponding to a situation when all active antineutrinos 
(or neutrinos) are converted into sterile neutrinos and to a value — 0.5 for 
the chemical potential Q. On the other hand it has been shown how, for tem- 

^ The authors of a recent paper obtain much lower values for the final lepton number 
flH]. This result has been obtained solving an approximated equation that, compared to 
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peratures T ~ 5 MeV, an effective Icinetic approach agrees almost perfectly 
with the quantum kinetic one for small mixing angles. It is also possible to 
extend the numerical study of QKE to large mixing angles |jl6[ and it has 
been found that for sin2 2^o ~ 10"^ and sin^ 2^o ~ 3 x 10-^(eVV|5m2|), in 
the quantum kinetic approach, at the critical temperature Tc, rapid oscilla- 
tions in the lepton number take place. This behaviour was first studied in 



1^ and later confirmed in [|T8| (even though in a smaller region of mixing 
parameters), but in these works the momentum dependence was not taken 
into account and the rapid oscillations are observed for a much larger region 
of mixing parameters. In [|16] is concluded that, outside the special region 
of mixing parameters where rapid oscillations are observed at the resonance, 
at temperatures T ~ 5 MeV, the effective kinetic approach provides a very 
good description not only for the evolution of the absolute value of lepton 
number but also in predicting its final sign. In this paper we will extend this 
analysis to the case when some tiny inhomogeneities in the baryon number 
(baryon isocurvature perturbations) are present, dealing with the region of 
mixing parameters where the effective kinetic approach can be used. This 
can be done including into the equation a term that describes the diffusion 
of neutrinos and we will show how this effect can induce the generation of 
lepton domains 0. 

the static approximation equation, contains a new term responsible, according to them, 
for the different values. We do not hide our scepticism toward these results, due to the 
robust coherent picture of previous ones, obtained both in a physically clear approximated 
picture ||, |l^, ^ and confirmed by numerical calculations using the exact QKE equations 
p3| , The authors of try to justify this situation claiming a not rigorous numerical 
procedure in the exact numerical calculations, without being aware of the conspiracy that 
would exist between a supposed numerical error in the solution of the exact QKE and the 
effect of a missing term in the static approximation equation. Moreover new numerical 
solutions of the exact QKE will be soon presented together with clear checks of their 
validity jl^ and an accurate description of the adopted numerical procedure. Although 
there is no final agreement at the present, we wish to thank A.D. Dolgov, S.H. Hansen, S. 
Pastor and D.V. Semikoz, for the kind availability to answer our questions. 

^ The possibility of a generation of lepton domains was first claimed in 0| due to a sign 
indetermination in the obtained equation for the lepton number evolution. This results 
in a generation of lepton domains with sign randomly determined. However in llo[| it was 
shown how the account of a correcting term produces a full sign determination of the 
solution (see Section 2). The idea of a "chaotic" generation of lepton domains has been 
recently re proposed in . This model assumes that the lepton number undergoes at the 
resonance very rapid and unstable changes of sign for any choice of mixing parameters and 
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The plan of the paper is the following. In section 2 we will show how the 
problem of sign of lepton number is fully determined in the effective kinetic 
approach, re analyzing and extending a procedure presented in |TD[. In sec- 



tion 3 we will extend the analysis to the case where small inhomogeneities in 
the baryon number are present, through the introduction of a diffusion term. 
We will illustrate how this term can induce, in some regions, an inversion in 
the sign of lepton number growth, bringing to the formation of lepton do- 
mains. In section 4 we describe qualitatively the evolution of lepton domains 
once they have been generated. In section 5 we conclude and discuss the 
possible applications of the new proposed mechanism. 



2 Final sign of lepton number in a homoge- 
neous background 

Within the Pauli-Boltzmann approach developed in |10|, the evolution of 
the lepton number carried by an active a-neutrino flavor that is mixed with 
a sterile neutrino flavor, is described by the following equation: 

"^^'^ [A{T,L)L-B{T,L)L,J (2) 



dt 



where we defined 



^ = ^ / dyy'fe,iT^, - f,,)(4 - 4) (3) 



"7 



6C(3)T=^ f , y 



,2 



B = ^'\^ / dy-^^—-{r^s + r^s) (4) 



7r%^ J ^ (1 + e2')2 



this would again result in a sign indetermination in different points of space. The analysis 
presented in including a full momentum dependence, excludes this possibility for 
almost all values of mixing parameters, except in a special region where the numerical 
calculations cannot be, at the moment, conclusive. In the mechanism we present here, the 
generation is not chaotic but perfectly determined by the spectral features of the baryon 
number inhomogeneities and moreover the horizon scale is not a limit to the size of lepton 
domains. 

^ We are actually neglecting a third term that arises only when a sterile neutrino asym- 
metry is produced, that means when the lepton number of active neutrinos is changing, 
considering that L^^ +L^^ =const. This term can give effects only when an initial a- lepton 
number much higher than L is assumed. Here we do not consider this situation and in 
this case this third term can be neglected. 
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The variable y in the integrals is the adimensional momentum p/T. The 
quantity L^^ is rigorously defined as the lepton number of the a-active neu- 
trino in the portion of comoving volume that contains a fixed number of 
photons A^*" at some initial temperature Tj„ ~ 150 MeV. Considering that 
the evolution of lepton number freezes at temperatures around T/ ~ 1 MeV, 
one can safely consider the number of photons in the element of comoving 
volume as a constant (neglecting muon annihilations) and write: 

Ny — Np — Uy 

We indicated with n the particle densities and with N = nR^ the numbers 
of particles in the comoving volume R^. We also introduced the (effective) 
total lepton number L defined as: 

L = U^+ + L,^ + L,^ ± = 2L,^ + L (6) 

where + holds for a = e and — for a = fi,T. The field L is the total charge 
number of non oscillating neutrinos plus a contribution from the baryon 
number carried by neutrons: we will refer to it as the background charge. It is 
constant while the oscillations occur, and must be considered as a parameter 
given by some earlier phase of baryo-leptogenesis. In this section we will 
assume that it is also strictly homogeneous, while in the next section we will 
study the effect of the presence of inhomogeneities. Both in L and in L we 
dropped an index a to simplify the notation. 

The quantities Tas, ^as are the production rates for sterile neutrinos and 
antineutrinos, f^g is the Fermi-Dirac distribution with zero chemical poten- 
tial and = {fua,s ^ fya,s)l'^feq the sum and difference distributions of 
neutrinos, relative to the Fermi-Dirac one. The active neutrinos distributions 
can be safely described by thermal equilibrium distributions because the pro- 
cess of generation of lepton number occurs for temperatures T 3> IMeV (the 
total collision rate is therefore much higher than the expansion rate) and for 
small mixing angles (the collision rates that refill the quantum states of active 
neutrinos are much higher than the sterile neutrinos production rates that 
deplete them). On the other hand the sterile neutrinos distributions must 
be described by two other rate equations that, together with the expressions 
for the production rates, can be found in 0, |ll| (we do not need them in 
the present context). 
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For a qualitative understanding of the evolution of lepton number, with 
a specific attention to its sign, it is useful to recast the equation (|]) in the 
form: 

^ = [2A(T, L) - B{T, L)] {U^ - L,,) (7) 
where we introduced the quantity: 

Leq = ^ (8) 

^ 2A-B ^ ^ 

While the term B is always positive, the term A changes sign at a critical 
temperature Tc, being negative for higher temperatures and positive for lower 
temperatures. When the temperature is far from the critical value (|T— ^ 
AT ~ (2 — 3) MeV), the B term acts as a correcting term [B ^ \A\) and 
the fixed point Lgg ~ L^^ = — (1/2)L. Moreover, until the lepton number 

is below a threshold value ~ 10~^(|5m^|/eV^)^ (see ^^), its non linear 
effect inside A and B can be neglected. In this situation the equation (|^) 
becomes extremely simple: 

^ = A{T) L (9) 

While for temperatures above Tc the quantity A is negative, the fixed point 
Leg is stable and thus the total lepton number is destroyed, for temperatures 
below Tc, A is positive, the fixed point is unstable and the total lepton number 
starts to grow. To answer the question toward which direction it grows, one 
has to take into account the term B and study its action in the vicinity of the 
critical temperature when A is small and B becomes important in driving 
the evolution of lepton number. Let us see in different stages what happens 
when temperature approaches the critical value from the above. When B is 
still small compared to A, expanding up to the first order in B/A, one gets: 

Leq = -\l + (10) 

During this period the term 2A — 5 is negative and the solution tracks the 
fixed point Lgg. An approximate expression for the quantity 5L = Ly^ — Lf.q 
can be derived neglecting a term d{5L)/dt in the Eq. (0): 

5L=—^ ^ (11) 

2A-B dt ^ ' 
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Until 2A — B ^ dLeq/dt, the solution tracks the fixed point. From the 



expression (|TOD it is clear that the term B drives the growth of the lepton 
number toward the same sign as the background charge L. At the critical 
temperature the term A vanishes and simply L^q = 0. In this moment the 
fixed point is changing very rapidly while the term 2A — B, that should 
force the solution L^^ to track L^g, is small and the solution starts to diverge 
from the fixed point. Immediately below the critical temperature, when 
2 A = B, the fixed point has a vertical asymptote and changes sign and at 
still lower temperatures it rapidly approaches zero again. This time however 
A is positive and the fixed point Leg does not attract the solution any more 
and this continues to grow toward the same direction transmitted by the 
"derailment" action of the term B at the critical temperature. 

The behaviour of the solution around the critical temperature is shown 
in fig. 1 for a particular choice of the mixing parameters. 



3 Lepton domains formation in presence of 
inhomogeneities 

In this section we will generalize the process of active-sterile neutrino os- 
cillations to the case when spatial inhomogeneities are present. From the 
observation of CMB we know that small inhomogeneities in the temperature 
field (adiabatic perturbations) were present in the early universe. However 
taking into account the presence of these perturbations does not change the 
basic results of the homogeneous scenario. The lepton number growth starts 
at slightly different times in different points, but the final result is unchanged: 
the presence in the end of a final lepton number with the same sign every- 
where in the space and with the same absolute value as in the homogeneous 
case. 

More interesting is to consider the possibility that small inhomogeneities 
are present in the background charge L (isocurvature perturbations). This 
quantity is the sum of the lepton number carried by the non oscillating neu- 
trinos and a term given by the presence of a baryon number. Whatever is the 
mechanism that created the inhomogeneities, it is reasonable to think that 
the inhomogeneities in the lepton numbers carried by the non oscillating neu- 
trinos become soon much smoother than those in the baryon number, due to 
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the much higher diffusion of neutrinos. More precisely we can consider the 
baryon number field expressed in the comoving coordinates as a constant, 
neglecting the neutron diffusion. This will make much simpler our following 
discussion, without altering the main results. 

A first trivial effect on the equation (Q) is simply that now L = L{x), 
where x are the comoving coordinates. At temperatures above this effect 
is enough to understand how the process is modified by the presence of the 
inhomogeneities. The Ua-lepton number will evolve in a way to destroy in 
any point the total lepton number L and a zero order solution, considering 
only the A-term, is given by L^^ = — (l/2)L(x). When temperature drops 
down, approaching the critical value, again the action of the correcting term 
B must be considered in order to understand toward which direction, positive 
or negative values, the total lepton number will grow. 

This time however a new correcting term must be considered in the equa- 
tion, a term that arises due to the neutrino diffusion, and the full equation 
becomes now: 



dL 



dt 



V 



A(T, L) [2L,„ + L(x)J - B(r, L) L,^ + -^^'L,. (12) 

where clearly L^^ = L^^{-x,t), V is the diffusion coefficient and R is the 
scale factor that appears because we are dealing with comoving coordinates 
(we normalize R in a way that Rq = 1), while the diffusion term must be 
calculated in physical lengths. 

The appearance of a diffusion term is something intuitively clear. For a 
more rigorous derivation one has simply to write the Liouville operator in 
the left hand side of the Boltzmann equation in the inhomogeneous case: 

d d c 

■^fuA^,P,t) = — /^,(x,p,t) - Hp-\/pfu^{x,p,t) + ■ Vx/i.„(x,p,t) 

(13) 

where H is the expansion rate R/R. In this way, integrating on the momenta 
the Boltzmann equation for the distribution difference of active neutrinos 
= fu^ — fp^ and dividing for the photon number density, one obtains: 



dL^ 
dt 



dL^ 
dt 



osc 



+ -^V.-(P V.L.J (14) 



where the first term in the r.h. side is the contribution from the oscillations 
into the sterile neutrinos (the r.h. side in the Eq. (^), while the second term 
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is the contribution due to the diffusion. The diffusion coefficient is defined 
by the expression: 

/(0y/.:(-.P.«)p--|v. / ^/.;(x,p.O (15) 

As we are considering the case of small inhomogeneities in the baryon num- 
ber, the diffusion coefficient can be safely considered homogeneous and one 
gets the diffusion term in the form written in the equation ([T^) . 

The order of magnitude of V is given by dint, where imt = ^ird is the 
interaction length and Tint = 3.15 /cqGI^T^ is the total collision rate, with 
ka - 1.27 (0.92) for a = e (/i, r) 0. 

If the diffusion term is negligible compared to the B term, then the ho- 
mogeneous scenario is practically unchanged. In the equation (|TD|) one has 
simply to consider that now L = L{x.) and, as we are considering small in- 
homogeneities so that the sign of L is spatially constant, the B term pushes 
the solution toward the same direction in all the points. 

On the other hand if we assume that the diffusion term is dominant, 
the situation can be much different. Let us assume for definiteness that 
L > 0. In the regions where the background charge is lower, the a-neutrino 
lepton number, before the critical temperature, would be higher (because 
L^^ ~ — 0.5L). In these regions the diffusion term pushes the lepton number 
to be depleted with the result that at the critical temperature one can have 
< —0.51/ and the growth starts toward a negative sign. The vice versa 
would happen in the regions where the background charge is higher. The 
final result is the creation of regions with different sign of lepton number 
with the same size as the scale of baryon number inhomogeneities. The sign 
of lepton number is in fact determined at any point by the curvature of the 
background charge field and thus the global properties of the baryon number 
inhomogeneities are somehow transmitted to the lepton domains. 

We have now to study the conditions for which this mechanism can work. 
The involved parameters are the mixing parameters, the amplitude of the 
inhomogeneities and the size of the inhomogeneities. We can assume, for 
simplicity, that inhomogeneities have only one size scale. Let us introduce 

An explicit calculation of the diffusion coefficient in the simpler case of " light" particles 
diffusing in a "heavy" particles medium can be found in pit]. 
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the field p(x) of tlie inliomogeneities, writing tlie background charge field as: 



L(x)=L[l + p(x)] (16) 

where L is the mean value, and where we assume that the field p is a pertur- 
bation (p -C 1 at any point). As we already said, while also non oscillating 
neutrinos give a non negligible contribution to the mean value, the baryon 
number gives the dominant contribution to the inhomogeneities. 

At temperatures higher than the solution will approach again the fixed 
point that at the zero order will be simply L^g{^) = — (1/2)1/ (x). The final 
sign of lepton number will be again determined by the correcting terms that 
this time include also the diffusion term. This can be generically estimated 
through the following expression (we drop the dependence on x in p and L'^g)'- 



V 



^^V%„^-r,.,p(^^j L%^-DL% (17) 

where we indicated with A the scale of the inhomogeneities and, given the 
generic physical length we mean with £(0) = (Ro/R) i, the comoving length 
normalized at the present. We also introduced the convenient quantity D. 
In this way the expression (0), valid for the fixed point in the homogeneous 
case, becomes: 

Le, = — (18) 

or at the first order in [B -\~ D)/A: 

= (19) 

This expression is now describing what we previously said in words: if the 
term D is negative (it means p < 0, the regions where the background charge 
is lower than the mean value if we assume it is positive) and its absolute 
value higher than B, then the lepton number growth is addressed toward the 
opposite sign than in the regions where the opposite condition holds. 

To determine the conditions for which an inversion of sign is possible, we 
have thus simply to compare the term B with the absolute value of D at 
temperatures around the critical temperature. An estimation of the order 
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of magnitude of B can be easily done considering that the integral on the 
momenta receives a dominant contribution only around the resonant value. 
This procedure has already been employed in []TT| (in that paper it was used to 



estimate the A-teim) and here we only give the result (valid for temperatures 



T ~ Tc, when the lepton number is still small and can be neglected): 

B ~ 10^ r,„i (20) 

where we defined s = sin26'o. A calculation of the comoving interaction 
length expressed in parsec gives the result: 

4«=. 100(70)pc(Mp^)' (21) 

for a = fi,T{e). Using this expression and imposing the condition \D\ > B 
for the generation of lepton domains, it is thus straightforward to derive the 
following condition on the (orders of magnitude) of the involved parameters: 



This analytical expression is confirmed by the numerical calculations. In 
figure 2 we show the evolution of the terms |y4|, B, calculated numerically for 
a particular choice of the mixing parameters and compared with the term \D\, 
calculated through the expression (|T^ with different values of the parameter 
\p\/{\^^^y. The result is in agreement with the one that could be derived 
using the eq. (|2^). 

This condition can also be expressed as a condition on the order of mag- 
nitude of the maximum size of lepton domains that can be generated 0: 



A(°) ~ lOKpc 













_10-5eV2_ 


vio-v 





< 



lOKpc (23) 



In this last expression we indicated the extreme possible values of the different 
quantities that yield the maximum value for A*-''-'. In particular the condition 

^ There is also a condition on the minimum size of lepton domains that can be gen- 
erated H: A(") ~ l^^liT ~ Tc) - 10-^pc{\dm'^\/eV'^yi. Below this scale neutrino are 
free streaming and they can destroy the lepton domains more rapidly than they can be 
generated, considering that L ~ Tint- 
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~ 10^^° must be imposed to have a not negligible final absolute value of 
lepton number. 

In next section we will use this result to sketch a scenario of how lepton 
domains with different sizes evolve and affect the BBN. This analysis will 
suggest a way to circumvent the upper limit found on the maximum size of 
a lepton domain. 

In conclusion of this section we notice that because of the presence of 
the term B there are always points where, even though p < 0, the condition 
B — \D\ < is not verified. Therefore we can say that the B term will 
always favour a dominance of regions where the generation of lepton number 
occurs with the same sign of L: we will refer to it as the dominant sign, 
compared to regions with inverted sign. There are however two much different 
extreme situations of dominance. In a first one B/\D\ <^ 1 everywhere, 
except in thin walls around the surfaces where p — {weak dominance). In 
this case the size of lepton domains coincides with the size of baryon number 
inhomogeneities and in a comoving sphere with a much larger radius than 
the size of inhomogeneities, the volume of regions with dominant sign lepton 
number is almost equal to that of regions with inverted sign lepton number. 
The topology of lepton domains is that one of a cubic lattice with each cube 
surrounded by first closest neighbours cubes with opposite sign and second 
closest neighbours cubes with the same sign. 

In a second extreme situation the condition B — \D\ < is verified only 
in small regions around the points where \p\ is maximum, with a size much 
smaller than that of the inhomogeneities {strong dominance). In this case 
one has a structure of lepton domains with inverted sign (3-dim) islands in 
a dominant sign background. 

4 Lepton domains evolution 

Lepton domains start to be formed at T = Tc. With the increase of lep- 
ton number the diffusion at the border of domains also increase. There is 
some interplay between the rate of generation (the A-term) and the rate of 
diffusion. 

Let us assume that the dominant sign is positive. In this case the diffusion 
can be described as a process that gradually fills at the border the regions 
with negative sign (the "holes"). At the same time, in the positive sign 
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regions, the lepton number that has been used to fill the holes, is restored by 
the oscillations through the growth term. In this situation it is clear that if 
the diffusion is able to cover the whole hole, the generation of lepton domains 
is only a transient regime without any practical effect. We can refer to this 
phase as cannibalization regime, in which the dominant sign regions enlarge 
at the expense of the inverted sign regions, until eventually their complete 
disappearance. 

However this process can last only until the generation of a lepton number 
can occur at the border of domains, where there is a change of sign and the 
lepton number is kept small by the diffusion. When the temperature drops 
down to T ~ 0.6 Tc or in any case down to 3 MeV, the growth starts to 
be inhibited. We have to calculate the diffusion length to know how large 
the lepton domains with inverted sign must be in order to survive to the 
cannibalization regime. This is given by the following expression that takes 
into account the Universe expansion: 

,(0)_ f'Vidt' 

The diffusion velocity is approximately given by Vd ~ \/T> jit — t(). Consid- 
ering that t ^ tf. and that t ~ it is an easy task to get the following 
expression for the comoving diffusion length: 

Ai°^ ~ 100(70) pc 

for a = fi,T (e). The generation of lepton number can only occur down 
to temperatures T ~ min(3 MeV, 0.6 Tc), therefore for ~ 5 MeV, using the 
expression (|l|), we get: 

Ar(T = 0.6Tc)~0.1pc(M^)"" (26) 

^ It corresponds to a situation when the resonant neutrinos have a momentum p ^ 10 T, 
in the tail of the distribution. For small values of L, neutrinos and antineutrinos are both 
resonant but if L > (< 0) and ii y > Upeak — 2.2, the number of resonant antineutrinos 
(neutrinos) is a little bit higher than the number of resonant neutrinos (antineutrinos): 
this explains why the fixed point L ~ is unstable and lepton number starts to grow (see 



^eV y I^MeV^ 
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For sizes below this value {very small scales) lepton domains with inverted 
sign are destroyed before they can produce any effect. For Tc ~ 5 MeV the 
diffusion length between the beginning of the generation of lepton number at 
Tc and its end at T ^ 3MeV tends to zero for |5m^| — > 10~^eV^, but in any 
case there is a lower limit on the size of lepton domains that can be generated 



2^ and in the end the same expression (^) can be approximately used. 

For larger scales the cannibalization cannot destroy completely the lepton 
domains before that the generation of lepton number at the border stops. 
When this happens, the lepton number that diffuses to fill the hole is not 
generated any more: the result is an effect of dilution of lepton number. The 
generation of lepton domains leaves a mark that modifies the homogeneous 
scenario. 

In this case, to understand the fate of lepton domains, it is useful to 
compare the scale of lepton domains with the horizon scale, given by the 
following expression: 

(f = lOOpc (M|X) (27) 



It coincides with the diffusion length at IMeV when neutrinos start to free 
stream. All lepton domains with small scales X^^^ ~ lOOpc, do not survive 
and one gets in the end a homogeneous lepton number field. However, as al- 
ready stated, this time the generation of lepton domains produced the effect 
to dilute the final value of lepton number. If the condition for the formation 
of domains (^) is satisfied only in the peaks of the inhomogeneities {strong 
dominance) , then the contribution from the regions with inverted sign is neg- 
ligible and the dilution effect too. Otherwise in the other extreme case, when 
the relation (^) is satisfied in almost all the space {weak dominance), then 
there would be an almost total reciprocal cancelation of regions with oppo- 
site sign, with only a small relic value of lepton number due to the presence 
of the 5-term. All intermediate values of lepton number are possible. 

Scales of inhomogeneities A^^-* ~ lOOpc can produce lepton domains able 
to survive until the freezing of neutron to proton ratio and therefore the 
presence of these scales would give rise to an inhomogeneous BBN scenario. 

Another important distinction is between scales able to produce primor- 
dial nuclear abundance inhomogeneities that survive until the present {large 
scales) or not {medium scales). In fact, even though nuclear abundances 
are produced in a inhomogeneous scenario of BBN, subsequent astrophysical 
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mixing mechanisms, such as shock waves induced by supernovae explosions, 
would be able to homogenize the products of BBN . These processes are 
effective on scales A*^°^ ~ lOOKpc. From this point of view the mechanism 
we proposed, in its simplest version, is unable to produce visible primordial 
abundance inhomogeneities. 

We can however circumvent this limit if we relax our initial assumption on 
the presence of only one characteristic scale in the field of the inhomogeneities 
p(x). We can in fact imagine a simple extension with two characteristic scales 
contemporarily present (for examples two Fourier modes): one small scale 
and one large scale. The amplitude of the first one could be modulated by the 
presence of the second. In this case one can have the simultaneous existence 
of regions where the small scales inhomogeneities are able to efficiently dilute 
the lepton number generated at the critical temperature down to negligible 
values, and regions where the amplitude of small inhomogeneities is depressed 
and there is no dilution effect. In this way one can produce islands where 
lepton number is present in a background with zero lepton number or vice 
versa. This time the size of the island regions has no intrinsic limit, but it is 
a feature of the model that produced the baryon inhomogeneities. 

More realistically one has to think that the inhomogeneities are described 
by a spectrum of lengths. The interesting property is that the spectral fea- 
tures of the isocurvature perturbations would be reflected in the lepton do- 
mains features and eventually in the BBN products. Another interesting 
aspect is that very large scale lepton domains would not be constrained 
by microwave background observations as models where large amplitude 
isocurvature perturbations are directly present in the baryon number 



This because the inhomogeneity in the active neutrino energy density would 
be compensated by an opposite inhomogeneity in the sterile neutrino en- 
ergy density. In fact during active-sterile neutrino oscillations one has that 
Ly^ + Ly^ =const. Even though sterile neutrinos free stream while active 
neutrinos are still diffusing, for very large scales (surely larger than the dif- 
fusion length of active neutrinos) this different behaviour cannot change the 
conservation of the sum of the asymmetries of active neutrinos and sterile 
neutrinos at each point. 

In the end of this section we want also to show that small amplitudes 
of inhomogeneities are sufficient for the mechanism to work. From the ex- 
pression (p2D , imposing the lower limit (|26|) on A^°^ (in order to have lepton 
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domains with inverted sign large enough not to be cannibahzed) , one gets 
the following condition on the mixing parameters: 

''^(^)'<10-V (28) 

that is satisfied for large regions of possible values of the mixing parameters, 
even for very small values of p. These regions are represented in figure 3 for 
three different values of p. It can be noticed that even for p = 10~^, there 
is a significant region where lepton domains can be generated and produce 
interesting effects. 



5 Conclusions 

We proposed a mechanism according to which active-sterile neutrino oscil- 
lations would amplify small baryon number inhomogeneities generating a 
structure of lepton domains. The assumption on the presence of baryon 
number inhomogeneities is a very reasonable one, because very small am- 
plitudes are sufficient to produce lepton domains large enough to give some 
effect. The structure of the lepton domains and its effect on the BBN would 
reflect the spectral features of the same perturbations that seeded the for- 
mation of the domains. The main attraction of this mechanism is that the 
size of these domains is not limited by the horizon scale. Moreover large 
scale lepton domains cannot be ruled out by current CMB observations (but 
maybe could be constrained or discovered in future experiments). 

The mechanism presents different applications. It provides a scenario for 
a non standard BBN (see for a review and see also the recent paper p6 



where the effects of inhomogeneous chemical potentials on BBN are consid- 
ered) . If the active neutrino is an electron neutrino the effects on the nuclear 



abundances can be remarkable [27, 28] and some of them could help to ex- 



plain the present observational picture, for example claimed inhomogeneities 
in the D/H abundances from measurements in quasar absorption systems 
[p9|] . From this point of view one possibility is that we live in a background 
of zero electron lepton number but some observed high redshift absorption 
systems, for example the one observed with an high abundance, could be 
included in an island of non zero lepton number Q 

^ The possibility that neutrino chemical potentials inhomogeneities can be responsible 
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We also mention that isocurvature perturbations can have interesting 



effects on large scale structure |32, 33 



Very generically we can say that the generation of lepton domains enlarge 
the possibility of the existence of some observational signature of active- 
sterile neutrino oscillations in the early universe. In the very fortunate case 
that active-sterile neutrino oscillations, with the right parameters for a lepton 
domain formation, do really occur in nature, then we would have a powerful 
probe for baryogenesis models that predict a spectrum of inhomogeneities. 
The simple non observation of a signature of the generation of lepton domains 
would put strong constraints. We conclude saying that the planned earth 
experiments will be able in next future to test the hypothesis of active-sterile 
neutrino oscillations and thus to rule out or support the proposed mechanism. 



Acknowledgements 

P. Di Bari acknowledges many interesting and stimulating discussions 
with R. Foot, P. Lipari, M. Lusignoli and R.R. Volkas and is grateful to them 
for many suggestions that made this manuscript clearer. He also thanks K. 
Jedamzik for a nice discussion during COSMO-99. 



References 

[1] R.Foot, M.J.Thomson and R.R. Volkas, Phys.Rev.BSS R5349 (1996). 

[2] K.Enqvist, K.Kainulainen and J.Maalampi, Nucl.Phys. B349 754 
(1991). 

[3] R.Barbieri and A.D. Dolgov, Nucl. Phy S.B349 743 (1991). 
[4] A.D. Dolgov. Sov.J.Nud.Phys. 33 700 (1981). 

[5] L. Wolfenstein. Phys.Rev. D17 2369 (1978). 

for the primordial deuterium abundance inhomogeneities able to explain the observations 
from high redshift quasar absorbers, has been proposed in | |30| . A recent analysis of 
the current measurements of primordial deuterium abundance in high redshift quasars 
absorbers can be found in |31|1. 



18 



[6] R.A.Harris and L. Stodolsky, Phys.Lett. B116 464 (1982). 
[7] L. Stodolsky, Phys.Rev. D36 2273 (1987). 

M.J.Thomson, Phys.Rev.A45 2243 (1992). 
[9] B.H.J.McKellar and M.J.Thomson, Phys.Rev.'D49 2710 (1994). 
[10] R.Foot and R.R.Volkas, Phys.Rev.'D55 5147 (1997); 
[11] P.Di Bari, RLipari, and M.Lusignoh, |hep-ph/9907548 . 



[12] Nicole F.Bell, Raymond R.Volkas, and Yvonne Y.Y.Wong, Phys.Rev. 
D59 113001 (1999). 

[13] R.Foot and R.R.Volkas, Phys.Rev.B56 6653 (1997); Erratum 
Phys.Rev.'D59 029901 (1999). 

[14] R. Foot, Astropart.Phys. 10 (1999). 

[15] A.D.Dolgov, S.H.Hansen, S.Pastor and D.V.Semikoz, |hep-ph/ 99 10444 

[16] P.Di Bari and R.Foot (in preparation). 

[17] X.Shi, Phys.Rev.TihA 2753 (1996). 



[18] K.Enqvist, K.Kainulainen and A.Sorri, [hep-ph / 9906452 



[19] X. Shi and G.M.FuUer |astro-ph/990"404l 



[20] K.Enqvist, K.Kainulainen and M.Thomson, Nucl.Phys.B373 498 
(1992). 

[21] L.D Landau, E.M. Lifsits, Physical Kinetic, (Pergamon, New York, 
1981). 

[22] R. Foot, private communication. 

[23] K.Jedamzik and G.M. Fuller, Ap.J. 452 33 (1995). 

[24] C.J.Copi, K.A.Olive and D.N.Schramm, Ap.J. 451 (1995), 

[25] K.Jedamzik, Proceedings of Neutrino Astrophysics, Ringberg, 1997. 

19 



[26] S.E.Whitmire and R. J.Scherrer, |astro-ph/ 99 10352 

[27] R.Foot, Ih^ph/ 99063 111 - 

[28] X.Shi, G.M. Fuller and K.Abazajian, Phys. Rev. D60 063002 (1999). 

[29] J.K.Webb et al., Nature 388 250 (1997). 

[30] A.D.Dolgov and B.E.J.Pagel, New Astronomy, in press. 

[31] D.Kirkman, D.Tytler, S. Buries, D. Lubin and J.M. O'Meara, ^stro- 

I ph/9907T2^ . 



[32] RJ.E. Peebles Nature 327 210 (1987). 

[33] R.Cen, ROstriker, and RJ.E. Peebles, Ap.J. 415 423 (1993). 



20 



Figure Captions 



Fig. 1 The behaviour of L^^ — L^^ (solid line) around the critical tem- 
perature when the sign is determined by the action of the correcting term B. 
The dotted line is the quantity Leg — L'^g, with Leg defined by the expression 
(H) and = — (1/2)L. We used L = 10~^°, while the values of mixing 
parameters are those indicated. 



Fig. 2 Comparison of the quantities A,B,D times \dt/dT\ = 1/HT (for 
numerical calculations it is more convenient to replace time with tempera- 
ture). The term A vanishes at the critical temperature T ~ 15MeV. The 
values of mixing parameters are those indicated. The diffusion term is plot- 
ted for different values of the parameter p/(A*-°'' pc^^)^. It can be noticed that 
the diffusion term |-D| is dominant on B, around the critical temperature, for 
p/(A^°^ pc"^)^ ~ 10~^, as it could be deduced by the expression (^2]). 



Fig. 3 Represantation of the regions where the condition (^) for the 
surviving of lepton domains to the canibalization regime is satisfied. The 
regions are those below the three solid lines and correspond to the three 
indicated values for the amplitude p of the baryon number inhomogeneities 
that seed the lepton domain formation. Above the dotted line the contri- 
bution of sterile neutrinos to the number of effective neutrinos during the 



BBN is higher than 0.8 for a = e (see pO[). Even for very small amplitudes 



(p = 10 ^) there is a significant region of parameters where the condition 
(EsD is satisfied. 
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